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ANALYTICAL  APPROXIMATE  SOLUTION  OF  COUPLED  WAVE 
EQUATIONS  WITH  A  NONLINEAR  STIFFNESS 


INTRODUCTION 


PURPOSE 


Equations  of  motion  can  be  derived  for  complex  dynamical  systems  by  using  calculus  of 
variations;  however,  because  of  the  nonlinear  complexities  of  these  equations,  closed-form 
solutions  are  usually  unattainable.  An  alternative  approach  is  to  interrogate  the  derived 
equations  by  numerical  methods.  Although  the  numerical  methods  do  allow  the  capture  of  the 
dynamic  response  of  the  system,  a  closed-form  solution  is  always  more  desirable.  This  report 
describes  the  Naval  Undersea  Warfare  Center  (NUWC)  Division  Newport’s  effort  to  develop  a 
method  to  derive  an  approximate  closed-form  solution  of  the  coupled  wave  equation. 


BACKGROUND 


In  the  early  1990s,  Adomian  developed  a  method  to  derive  analytical  approximate 
solutions  to  nonlinear  functional  equations.1  This  method  is  referred  to  as  the  “Adomian 
decomposition  method”  (ADM).  The  solution  to  the  given  nonlinear  functional  equation  can  be 
approximated  by  an  infinite  series  solution  of  the  linear  and  nonlinear  terms,  provided  the 
nonlinear  terms  are  represented  by  a  sum  of  series  of  Adomian  polynomials.  ’  ’  ADM  has  been 
successfully  applied  to  various  types  of  ordinary,5’ 6' 7  partial,8’ 9’ 10’  n’ 12  and  delay  differential 
equations  to  develop  closed-form  approximate  solutions. 


SCOPE 

In  this  report,  the  derivation  of  an  approximate  closed-form  solution  of  the  coupled  wave 
equation  is  accomplished  by  using  the  ADM.  The  coupling  is  realized  by  both  the  nonlinear 
softening  and  the  nonlinear  hardening  springs  and  a  linear  spring  located  at  the  center  of  the 
crossing  strings.  The  dynamic  response  of  the  closed-form  solutions  is  compared  for  each  spring 
type  under  the  same  prescribed  initial  condition. 

The  remainder  of  this  report  explores  the  Adomian  decomposition  method  and  provides  a 
brief  derivation  of  the  coupled  wave  equations;  derives  approximate  analytical  solutions  for  the 
coupled  wave  equation  using  the  Adomian  decomposition  method;  and  investigates  the  closed- 
form  solutions  for  the  three  different  types  of  springs  under  the  same  prescribed  initial  condition. 
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ADOMIAN  DECOMPOSITION  METHOD 


GENERAL 

Consider  equation  (1), 

Fu  =  g(t),  (1) 

where  F  is  a  nonlinear  operator,  and,  when  it  is  expanded  into  linear  and  nonlinear  terms, 
results  in, 

Lu  +  Ru  +  Nu  =  g([).  (2) 

The  linear  term  is  represented  by  Lu,  L  is  invertible  and  is  chosen  to  be  the  highest  order 
derivative,  R  is  the  remainder  of  the  linear  terms,  g(t)  is  the  nonhomogeneous  right  side,  and 

the  nonlinear  terms  are  represented  by  Nu.  If  L  is  chosen  to  be  a  derivative  of  order  n ,  then  L  1 
will  be  an  n  -fold  integral.  For  example,  if 

d1  1  ' 

L - -  [•]  then  Ll  =  J  J  [\itdt, 

dt  oo 

and 

r1L  =  u(0)+tu'(0). 

Solving  equation  (2)  for  Lu  results  in 

Lu  =  -Ru  -  Nu  +  g(t).  (3) 

Adomian  postulated  that  the  solution  u  can  be  approximated  by  an  infinite  sum  of  series 

oo 

«  =  (4) 

n= 0 

and  the  nonlinear  term  Nu  can  be  approximated  by  a  sum  of  Adomian  polynomials 

oo 

Nu  =  ^4>0,Ui,...,w„) ,  (5) 

n= 0 

where  the  Adomian  polynomials  of  u(),u],...,un  can  be  calculated  by  the  formula 
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An(u0,u1,...,un)  = 


1  d" 


n !  dAn 


N 


C  ®  A 

v  1=0  y 


(6) 


A=0 


where  A  is  a  dummy  variable.  The  first  three  Adomian  polynomials  are  shown  for  convenience, 


A()  =  N(u0) 

4  =u1N'(u0)  (7) 

ir 

A2=-±N'(u0)  +  u2N(u0), 


where  (•)'  denotes  the  first  derivative  with  respect  to  u  .  Substituting  equation  (4)  and  equation 
(5)  into  equation  (3)  results  in 


oo  oo  oo 

n= 0  n= 0  n= 0 


(8) 


and  operating  on  equation  (8)  with  its  inverse  L  1  yields 

oo  oo  oo 

X  un  =  W(°)  +  tU '(°)  “  L'RYaUn  ~  L  '  X  An  (M0  >  Ul  ’ - 
«=0  «=0  n=0 

Finally,  equation  (9)  can  be  rewritten  in  a  recursive  relation  as  in  equation  (10), 

u0  =  u{  0)  +  tu\ 0)  +  Llg(t), 
u„  =  -L'Run_]-L'An_],  n>  1, 


(9) 


(10) 


where  u(]  is  obtained  from  the  prescribed  initial  and  boundary  conditions.  The  A-term 
approximation  can  be  used  as  a  series  solution, 

k-l 

k  =  u~  Yjui  ’  (H) 

/=0 

where  the  exact  solution  is  given  by 

u  =  lim  un.  (12) 

n—>  co 

If  more  accuracy  is  desirable,  then  more  terms  should  be  included  in  the  series  solution. 
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WAVE  EQUATION  DERIVATION 


Using  the  Euler-Lagrange  equations  allows  the  coupled  wave  equations  to  be  derived 
where  the  coupling  is  realized  through  a  nonlinear  softening  spring.  Consider  the  two  crossing 
strings  depicted  in  figure  1,  where  the  strings  he  in  the  xiX2-plane.  In  this  report,  the  strings  are 
considered  to  be  orthogonal  to  each  other  throughout  the  simulation.  Shearing  or  twisting  effects 
of  the  strings  are  not  considered.  The  string  in  the  x,  -direction  has  a  vertical  displacement 

u(xl,x2,t) ,  and  the  string  in  the  x2  -plane  has  a  vertical  displacement  v(xx,x2,i) . 


Figure  1.  Two  Crossing  Strings  Where  the  Coupling  Is  Realized 
with  a  Nonlinear  Softening  Spring 


The  Lagrangian  can  be  defined  as 

L(x1,x2,x1,x2,t)  =  T-Vf-Vs, 

where  T  is  the  kinetic  energy,  Vf  is  the  potential  energy  due  to  the  work  required  to  displace  a 
differential  element  of  the  string,  and  Vs  is  the  potential  energy  of  the  nonlinear  spring  where  a 

nominal  spring  length  is  used.  After  substituting  in  the  terms,  the  Lagrangian  can  be  expressed 
as 


L  =  —p 
2 


(d  V  (8  > 

—  w(Xj,X2,t)  +  — v(Xj,X2,t) 


2  X 


E_ 

'~2 


dt 


dt 


du{xx,x2,t)  du(xl,x2,t) 


V  dxi 


-  +  - 


ox 


+ 


2  J 


dv(x1  ,x2,t)  8  v(x,,x2,Q 


\2  V 


fix, 


-  +  - 


fix. 


2  J 


(13) 


( v(Xj ,  x2 ,  t)  -  m(Xj ,  x2 ,  t)\  -  ( v(Xj ,  x2 ,  t)  -  u(xl ,  x2, 0)4 
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where  p  is  the  mass  density  of  the  string,  E  is  the  Young’s  modulus  of  the  string,  and  K L  and 
Knl  are  the  linear  and  nonlinear  stiffness  coefficients,  respectively.  In  this  report,  p  and  E  are 

kept  the  same  in  both  strings.  In  its  normal  form,  the  Euler-Lagrangian  differential  equation  is 
expressed  as 


dL_d_fdL' 
dq  dt  ^  dq , 


(14) 


where  q  is  the  generalized  coordinate  and  q  is  the  time  derivative.  The  system  of  interest  here, 
however,  is  continuous,  so  the  generalized  coordinates  are  u(x],x2,t)  and v(x1,x2,t)  ;  thus,  the 
Euler-Lagrangian  differential  equations  for  u(x] ,x2,t)  and  v(x1  ,x2,t)  are 


dL 


d 


8u(x1,x2,t)  dt 


dL 


du(xvx2,t) 


dL 


d 


dv(x j ,  x2 ,  t)  dt 


dt 


dL 


dv{xvx2,t) 


d 


dt 


_d_ 

dx. 


_d_ 

dx , 


dL 


^  du(x  X,x2,t)^ 
V  dxi  J 


d 


dL 

^  dv(x j ,x2,t)'[ 


dx, 


i  J 


dx \ 


_d_ 

fix. 


dL 


^  du(xl,x2,t)'^ 


dx. 


2  J 


=  0 


dL 


^  dv(xx,x2,t 


V1 9  2 

dx. 


2  J 


(15) 


=  0. 


Evaluating  the  terms  yields  the  coupled  wave  equations 

d2u(xl,x2,t) 


dt2 


=  c2V2u(xt  ,x2,t)  +  KL  ( v(x,  ,x2,t)~ u(x j , x2 , t))  - Knl  ( v(xx ,x2,t)~ u(xt , x2 , t)) 


d  v(x,,x2,0  =  c2V2v^ ^t)  +  KL( u(x, , x, , t) - v(x, , x2, 0) - Knl  (u(xx ,x2,t)~ v(x, , x2 , t)f  , 
dt 


(16) 


where 


2  IE 


2  IE 


ci  =J~  and  c2  =■!—■ 


p 


p 
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ADM  APPLIED  TO  THE  COUPLED  WAVE  EQUATION 


According  to  the  ADM,  the  solution  to  u  -  u(xl,x2,t)  andv  =  v(xx,x2,t)  can  be 
represented  as 

co 

u(xl,x2,t)  =  YJun(xl,x2,t) 

and  "=0  (17) 

oo 

v(xvx2,t)  =  YJvn(xl,x2,t). 

n- 0 


Because  the  summation  of  un  and  vn  is  the  desired  solution  and  converges  rapidly,  the  k  -tenns 
of  the  summation  can  be  calculated.  The  resulting  approximate  analytical  solution,  therefore,  is 


*-i 

</>k  =XM»(X1’X2>0 

and  "=°  ( 1 8) 

%  =ZV»(X 1>X2»0- 

n= 0 


Because  the  nonlinear  terms  are  functions  of  two  variables,  one  cannot  continue  with  the 
Adomian  polynomials  described  in  equation  (6).  The  Adomian  polynomials  for  the  two 
displacement  variables  u0,ux,..., un , v0 , v, , . . . , vn  are  described  by 


Bn(u0, n„,v0,v1,...,v)i)  = 


1  d" 


nldr 


N 


V  i=0 


i=0 


A=0 


where  the  first  three  terms  are  shown  for  convenience: 


(19) 


B0=N(u0,v0) 

B^Ui—Niu^v^  +  Vy—Niu^Vo)  (20) 

ou  ov 


2  cs2 

2  2!  du 


v2  d2 

N(U0 ,  v0 )  +  ~7  N(u0  ,v0)  +  ulv] 


2!  dv~ 


^2  ^  q 

1  K 1  *(«0  ’ V0  )  +  U2  ~  N(U0  ’Vo)  +  V2^~  N(U0  ,  VQ  ). 

OUOV  (II  ov 


For  the  remainder  of  this  report,  u  -  u(xvx2,t)  andv  =  v{xvx2,t) .  From  the  coupled 
wave  equations  described  in  equation  (16),  the  highest-order  derivative  is  second  order,  so  the 

d2  ‘  1 

linear  operator  becomes  L  =  — r  and  L  1  =  j*  J  (-)dtdt .  In  the  u  equation,  the  nonlinear  tenn  is 


0  0 


N(u,  v)  =  Knl(v  -  u )3  and  N(u,v )  =  KNL  (n  -  v)3  for  the  v  equation.  The  first  three  terms  of  the 
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Adomian  polynomial  series  expansion  are  shown  for  convenience  where  the  superscripts  u,  v 
denote  the  polynomials  for  the  u,v  equations,  respectively: 


B‘‘=KNL(v0-u0f 

B['  =3KNL[(v0-u0)\vl-u1)~] 


=  Knl 


b; = knl 


3(Vj  -  i/j  f  (■ v0  -  u0 )  +  - (v0  -  u0  f  ( 2v2  -2 u2) 

3  1 

(Vj  -  w,)  +  3(v0  -  w0)(v1  -  Mj)( 2v2  -  2 u2)  +  -(v0  -  w0)  (6v3  -  6w3) 


(21) 


DV  _  n« 

5;  =  -5; 

DV  _  n« 

DV  _  n« 

z>3  —  —n 3 . 


The  /7  th  recursive  equation  for  both  w  and  v  are  found  by  plugging  equation  (16)  into  equation 
(10),  which  yields: 


u„  = 


1  1  1  1 

JJ[AV  i+KL{vn-l-un_l)ytdt-\\Bll(u0,...,un_l,v0,...,vn_l)dtdt 


0  0 
t  t 


0  0 
t  t 


(22) 


vn  =\\\cl^2vn.l+KL{un_l-vn_l)ytdt-\\Bl_l{u(),...,un_l,v0,...,vn_l)dtdt. 


Because  the  solution  is  a  sum  of  iterates,  it  must  be  transformed  into  a  closed-form 
solution.  This  transformation  is  achieved  by  taking  a  Laplace  transfonn  of  both  (f)k  and  cpk  with 
respect  to  t .  A  Pade  approximation  of  order  [2,  2]  was  used  to  approximate  the  coefficients  of 
the  Laplace  transformed  polynomials  into  a  rational  function.14  An  inverse  Laplace  transform 
was  used  to  convert  back  to  the  time  domain  to  yield  a  closed-fonn  solution. 
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SIMULATION  RESULTS 


In  this  section,  the  approximate  solutions  are  computed  for  three  types  of  springs  with  the 
same  initial  condition  with  the  following  system  parameters  over  the  domain  0 .  Both  strings 
experience  an  initial  displacement  ( u0  and  v0 )  as  shown  in  equation  (23)  with  zero  initial 

velocity  ( ux  and  vx ): 

u0  =  sin(7rxl)sin(;rx2),  v0  =  -^sin(;rx1)sin(;rx2),  m,  =  v,  =  0 
c,2=c22=  1,  g(t)=0,  0  =  [0,1] x [0,1]. 

The  following  spring  constants  were  chosen: 

(1)  nonlinear  softening  spring  KL  =  1,  K  v/  =  -3 

(2)  nonlinear  hardening  spring  K,  =  1,  K  V7  =  3  (24) 

(3)  linear  spring  KL  =  1,  K  Y/  =  0. 

In  the  first  case,  the  steps  used  toward  deriving  the  closed-fonn  solutions  are  presented;  however, 
in  the  other  two  cases,  only  the  responses  are  provided. 
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CASE  1 


In  this  case,  a  nonlinear  softening  spring  is  used.  Only  the  first  three  iterations  are  shown 
due  to  the  number  of  tenns;  four  iterations,  however,  were  used  to  develop  the  closed-fonn 
solution 


u0  =  sinf/rx,  )sin(;rx2) 
v0  =  -^-sini/TA',  )sin(;rx2 ) 


Wj  = - sm{nxx )  sin(;rx2 )  ( 1 6  n1  +93-81  cos2  (;rx2 )  -  8 1  cos2  {nx  \ )  +  8 1  cos2  ( nxx )  cos2  (;rx2 )) 


Vi  =r 


u2  = - sin(^x1)sin(^x2) 


192 


16 

sin(^Xj )  sin(;rx2 )  (8;r2  +  93  -  8 1  cos2  (;rx2 )  -  8 1  cos2  {nx \ )  +  8 1  cos2  {nx  \ )  cos2  (;rx2 )) 

22785  +  32;r4  -  3888^2(cos2(^x1)  +  cos2(^x2)) 

+8 1 8 1 0  cos2  {nxx )  cos2  (;rx2 )  +  3450;r2 
-42444(cos2(;rx1)  +  cos2(;rx2))  +  4374;r  cos2(;rx1)cos2(;rx2) 
-39366(cos4  (;rx2 )  cos2  {nxx )  +  cos4  {nx  x )  cos2  (;rx2 )) 
v+19683(cos4(^Xj)  +  cos4(^x2)  +  cos4(^Xj)cos4(^x2)) 

22785  +  32;r4  - 3888;r2(cos2(;rx1)  +  cos2(;rx2)) 

+8 1 8 1 0  cos2  {nx  \ )  cos2  (;rx2 )  +  3450;r2 
-42444(cos2(;rx1)  +  cos2(^x2))  +  4374;r  cos2(;rx1)cos2(;rx2) 
-393  66(cos4  (;rx2 )  cos2  {nx  x )  +  cos4  {nx  x )  cos2  (;rx2 )) 
+19683(cos4(;rx1)  +  cos4(;rx2)  +  cos4(;rx1)cos4(;rx2)) 


v2  = - sinf/rx,  )sin(7rx2) 


192 


(25) 


After  four  iterations,  the  approximate  solution  for  u  is  <j)4{xvx2,t)  =  'y'u,,  ;  for  v,  the 

n= 0 
3 

approximate  solution  is (pA (Xj ,x2,t)  =  ^ un  .  The  Laplace  transform  of  <f)A (x, ,x2,/)  and 


n= 0 

^4(xj,x2,t)  yield  </>4{xx,x2,s)  and<p4(xj,x2,s)  .  The  Pade  approximate  of  order  [2,  2]  for  both 
<fiA{xl,x2,s)  and  cpA{xvx2,s)  is 


9 


,  ,  .  81sin(^x)sin(^x,) 

<pd(x,,x,,s)  =  —7 - 1 - - -  N 

4  1  2  f  6561  /  2.  ,  2,  .  ,  2.  ^  ...  2  7533  ^ 

I  cos  (^"Xjjcos  (;rx2)-cos  (/rx,)-cos  {nx2) )  + 1  62tt  + 


81  + 


8 


8 


(pA{xl,x2,s)  =  ~- 


9sin(/rxl  )sin(/z\x2) 


(26) 


81  + 


^ cos2  {nx j )  cos2  (;rx2 )  -  cos2  {nxx )  -  cos2  (;rx2 ))  + 1 62 n1  + 


A ' 


The  final  closed-form  solution  is: 


u(xl,x2,t )  =  sin(/rx1)sin(/rx2)cosh^-|j-^-186-32/r2  +  162cos2(^x1)  -162cos2(^x2)(cos2(^x1)  - 1) 
v(Xj,x2,0  =  -  ^  sinl/TX, )  sin(/rx2 )  cosh  f  yj— 93  -  8/r  +  8 1  cos2  (xxl )  -  8 1  cos2  (/rx2  )(cos2  (/rx1 )  - 1)  j . 


\ 

(27) 


The  dynamic  response  to  the  prescribed  initial  condition  is  shown  in  figures  2  and  3. 
Figure  2  depicts  the  response  of  v(x,,x2,t)  (blue  line)  along  the  x2 -direction  and  u(xx,x2,t)  (red 
line)  along  the  x1  -direction  at  time  snaps  t  =  0.0,  0.4,  0.8,  1.2,  1.6,  and  2.0.  Figure  3  shows  a 
projection  of  these  two  responses  for  v(x,,x2,t)  (blue  line)  and  «(x,,x2,/)  (red  line)  along  the 
same  coordinate  axis  to  compare  the  oscillations  between  the  two  strings.  Initially,  the  spring  is 
compressed  with  the  string  along  the  x,  -direction  having  an  initial  displacement  with  a  greater 

magnitude  than  the  string  along  the  x2  -direction.  The  effect  of  the  nonlinear  softening  spring 

becomes  evident  as  time  progresses.  The  center  location  of  both  strings,  where  the  spring  is 
attached,  oscillates  at  a  higher  frequency  than  the  points  closer  to  the  fixed  ends.  In  essence,  this 
oscillation  creates  a  rippling  effect  from  the  center  of  the  string  to  the  boundaries.  If  time  were 
to  progress  further,  the  string  would  experience  high-frequency  oscillations  throughout  its  length. 
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t  =  0.00  t  =  0.40  t  =  0.80 


Figure  2.  Response  of  v(x1,x2,t)  (blue  line)  Along  the  x2-Direction  and  u(xl,x2,t) 

(red  line)  Along  the  x,  -Direction  for  Case  1  at  Time  Snaps  t  =0.0,  0.4,  0.8, 1.2, 1.6,  and  2.0 
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t  =  0.00 


t  =  0.40 


t=  0.80 


t  =  1 .20 


t  =  1 .60 


t  =  2.00 


Figure  3.  Projection  View  of  the  Response  of  v(xl,x2,t)  (blue  line)  and  u(xt,x2,t)  (red  line) 
Along  the  Same  Axis  for  Case  1  at  Time  Snaps  t  =0.0,  0.4,  0.8, 1.2, 1.6,  and  2.0 


12 


CASE  2 


In  this  case,  a  nonlinear  hardening  spring  is  used.  The  final  closed-form  solution  is: 


u(xl,x2,t)  =  sin(^x1)sin(^-x2)cosh^yl38  -32^-"  -  162cos  (ttXj)  +  162 cos' (;rx2)(cos  {nx j)- 1) 
v(Xj ,x2,t)  =  -—sin(7rxl)sin(7rx2)cosh^tyj69-87r2  - 81cos2(^Xj)  +  81cos2(^x2)(cos2(^x,) - 1)  j. 


(28) 


The  dynamic  responses  depicted  in  figures  4  and  5  are  similar  to  those  for  case  1  dynamics  in  figures  2 
and  3.  This  similarity  can  be  observed  in  the  closed-form  solutions  in  equations  (28)  and  (27),  which 
differ  only  by  constants.  The  global  oscillations  remain  the  same,  but  the  higher  frequency  oscillations 
are  different  from  what  was  observed  in  the  nonlinear  softening  spring.  Figure  4  depicts  the  response 
of  v(xj,x2,t)  (blue  line)  along  the  x2 -direction  and  the  response  of  u(x1,x2,t)  (red  line)  along  the  Xj- 
direction  at  time  snaps  t  =  0.0,  0.4,  0.8,  1.2,  1.6,  and  2.0;  figure  5  shows  a  projection  of  these  two 
responses  for  v(xl,x2,t)  (blue  line)  and  i<(x!,x2,f)  (red  line)  along  the  same  coordinate  axis.  The 
effects  of  the  nonlinear  hardening  characteristics  of  the  spring  are  apparent  in  the  response  of  the 
center  location  of  the  strings  where  the  spring  is  located.  Unlike  the  nonlinear  softening  case  (case  1), 
both  strings  in  case  2  have  their  inflection  points  at  different  times  in  the  simulation. 


t  =  0.00  t  =  0.40  t  =  0.80 


Figure  4.  Response  of  v(xl,x2,t)  (blue  line)  Along  the  x2 -Direction  and  u{xvx2,t)  (red  line) 
Along  the  Xj  -Direction  for  Case  2  at  Time  Snaps  t  =0.0,  0.4,  0.8, 1.2, 1.6,  and  2.0 
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t  =  0.00 


t  =  0.40 


t  =  0.80 


t  =  1.20 


t  =  1.60 


t  =  2.00 


Figure  5.  Projection  View  of  the  Response  of  v(x1,x2,t)  (blue  line)  and  u(xl,x2,t )  (red  line) 
Along  the  Same  Axis  for  Case  2  at  Time  Snaps  t  =0.0,  0.4,  0.8, 1.2, 1.6,  and  2.0 
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CASE  3 


In  this  case,  a  linear  spring  is  used.  The  final  closed-fonn  solution  is: 


v(x  \,x2,t)  =  - — sin^Xj )  sin(;rx2 )  cosh  [t'-j-ln:1  -  3  j . 


(29) 


The  dynamic  response  for  case  3  is  depicted  in  figures  6  and  7.  Because  the  coupling 
spring  connecting  the  two  strings  is  linear  and  the  initial  condition  is  composed  of  just  the  first 
mode  of  oscillation  for  a  string  with  fixed  ends,  both  strings  oscillate  only  in  that  first  mode.  The 
figures  show  that  there  are  no  higher  frequency  oscillations  in  the  string,  which  can  further  be 
seen  in  the  close-form  solution  in  equation  (29). 


t  =  o.oo 


t  =  0.40 


t  =  0.80 


Figure  6.  Response  of  v(x, , x2 ,  t)  (blue  line)  Along  the  x2  -Direction  and  u(xvx2,t)  (red  line) 
Along  the  Xj  -Direction  for  Case  3  at  Time  Snaps  t  =0.0,  0.4,  0.8, 1.2, 1.6,  and  2.0 
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t  =  0.00 


t=  0.40 


t  =  0.80 


> 

=r 


t  =  1 .20 


t  =  1 .60 


> 

3" 


0.5 

X1’X2 


Figure  7.  Projection  View  of  the  Response  of  v(xx,x2,t)  (blue  line)  and  u(x],x2,t)  (red  line) 
Along  the  Same  Axis  for  Case  3  at  Time  Snaps  t=  0.0,  0.4,  0.8, 1.2, 1.6,  and  2.0 
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CONCLUSIONS 


The  Adomian  decomposition  method  was  used  to  develop  an  approximate  analytical 
solution  for  coupled  wave  equations.  The  coupling  was  realized  by  nonlinear  softening  and 
nonlinear  hardening  springs  as  well  as  for  a  linear  spring.  The  equations  for  each  type  of  spring 
were  derived  from  the  Euler-Lagrange  equations.  The  effects  of  the  nonlinear  spring  behaviors 
were  evident  in  both  the  dynamic  response  and  closed-form  equations — very  unlike  the  behavior 
observed  for  the  linear  spring  case.  In  both  the  nonlinear  softening  and  nonlinear  hardening 
springs,  higher  frequency  oscillations  were  observed  at  the  spring  location;  however,  when  a 
linear  spring  was  used,  the  strings  exhibited  no  higher  frequency  oscillations  and  oscillated  only 
in  the  mode  by  which  they  were  excited. 
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